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The amplitude mode is a ubiquitous collective excitation in condensed-matter systems with broken 
continuous symmetry. It is expected in antiferromagnets, short coherence length superconductors, 
charge density waves, and lattice Bose condensates. Its detection is a valuable test of the correspond- 
ing field theory, and its mass gap measures the proximity to a quantum critical point. However, 
since the amplitude mode can decay into low-energy Goldstone modes, its experimental visibility 
has been questioned. Here we show that the visibility depends on the symmetry of the measured 
susceptibility. The longitudinal susceptibility diverges at low frequency as Xv<? ~ ioj^ 1 (d — 2) or 
log(l/|u;|) (d — 3), which can completely obscure the amplitude peak. In contrast, the scalar sus- 
ceptibility is suppressed by four extra powers of frequency, exposing the amplitude peak throughout 
the ordered phase. We discuss experimental setups for measuring the scalar susceptibility. The 
conductivity of the 0(2) theory (relativistic superfluid) is a scalar response and therefore exhibits 
suppressed absorption below the Higgs mass threshold, a ~ w 2 . In layered, short coherence 
length superconductors, (relevant e.g. to cuprates) this threshold is raised by the interlayer plasma 
frequency. 



I. INTRODUCTION 

A fundamental consequence of spontaneous (continu- 
ous) symmetry breaking (SSB) of ./V-component order pa- 
rameters is the emergence of collective order parameter 
oscillations: Goldstone modes and a massive amplitude 
(Higgs) mode. Examples in condensed matter systems 
are plentiful: for example, Heisenberg and XY spin sys- 
tems, superconductors, cold atom condensates in optical 
lattices and incommensurate charge density waves P In 
particle theory, the Higgs bosorP is modeled by the am- 
plitude mode of a gauged bosonic condensate. 

While long-wavelength Goldstone modes are sharp ex- 
citations in the broken symmetry phase, the amplitude 
mode is long-lived only at the classical (weak coupling) 
level. Quantum corrections allow for i ts d ecay into pairs 
of Goldstone modes. Previous author^El have therefore 
questioned the experimental visibility of the amplitude 
mode. Indeed, in ordered phases of two and three dimen- 
sions, the longitudinal susceptibility exhibits an infrared 
singularity, which can broaden the amplitude mode peak 
into an undetectable shoulder, as in Fig. [T] 

Is the amplitude mode therefore overdamped? This 
paper shows that, surprisingly, it is not. Rather, the 
infrared singularity is a property of the chosen suscep- 
tibility. The amplitude mode produces, in fact, a pro- 
nounced peak in the scalar (rotationally invariant) sus- 
ceptibility, as depicted in Fig. m Similarly, for the 0(2) 
model (relativistic supernuids), the mode produces a 
well-defined pseudogap in the optical conductivity, as 
previously noted in Kubo formula computations!^ Here 
we find at two-loop order an additional weak absorption 
tail below the Higgs mass threshold, a ~ w 2d+1 . 

The two types of susceptibility are associated with dif- 
ferent ways of parametrizing the fluctuations of the order 
parameter 4>. Within the ordered phase, the order pa- 
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FIG. 1. Visibility of the amplitude mode. The longitu- 
dinal and scalar susceptibilities of the two-dimensional O(N) 
model at zero temperature are plotted in the large N limit, 
m is the renormalized amplitude mode mass. The coupling 
constant is g — 0.84 g%° , which is within the broken sym- 
metry phase. Note that dissipation into Goldstone modes has 
very different effects on the low- frequency behavior of the two 
susceptibilities and on the visibility of the amplitude mode 
around uj — m. 



rameter gets an expectation value |(4>)| = 3>0j an d fluc- 
tuations in its value can be written as 

* = ($0 + 0",7r) 

( r~ \ . (!) 
= $ + VNp) h. 

In the first parametrization, the fields a and tv describe, 
in turn, the longitudinal and transverse fluctuations rela- 
tive to the ordering direction. In the second parametriza- 
tion, p is the fluctuation in the magnitude of the order 
parameter, while n is the ordering direction. We find 
that the scalar susceptibility, associated with p, can dis- 
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play a sharp peak at the Higgs mass, even in cases where 
the longitudinal susceptibility, associated with a, does 
not. 

The difference between longitudinal and amplitude 
fluctuations can be understood heuristically by consid- 
ering oscillations of a particle near the minimum of the 
Mexican hat potential shown in Fig. [2] The longitudi- 
nal component (in the broken symmetry direction) loses 
coherence rapidly as the particle oscillates and also me- 
anders around the rim. In contrast, amplitude oscil- 
lations, that is, fluctuations in the radial distance p, 
are much longer lived since their interaction with long- 
wavelength Goldstone fluctuations is suppressed by two 
derivatives. Thus, the infrared singularity in the longi- 
tudinal response is due to "contamination" of Goldstone 
modes in the response function and not due to overdamp- 
ing of the amplitude mode it self This explanation is 
made precise by calculating the scalar susceptibility of 
the relativistic O(N) field theory in the broken symme- 
try phase. 

This paper is organized as follows: In Sec.|TlJ we intro- 
duce the O(N) field theory on which we base our anal- 
ysis, define longitudinal and scalar susceptibilities of the 
order parameter, and show how to define conductivity in 



O(N) models. In Sec. Ill we compute in the weak cou- 



pling limit the susceptibilities to one-loop order and the 
optical conductivity up to two-loop order. We show that 
a cancellation of self-energy and vertex diagrams sup- 
press the sub-gap spectra of the scalar susceptibility and 
the 0(2) conductivity by four powers of frequency rela- 
tive to non symmetric response functions. In Sec. |IV| we 
compute the O(N) susceptibilities in the large N limit 
and find a suppression by the same four powers of fre- 
quency in the scalar susceptibility relative to the longi- 
tudinal susceptibility. We conclude in Sec. [V] by propos- 
ing various scalar probes for the amplitude mode in cold 
atoms, magnetic systems, and layered superconductors 
with long-range Coulomb interactions. In the cuprates, 
for example, we suggest that the mid-infrared spectral 
weight in the optical conductivity might be related to 
the amplitude mode of tightly bound Cooper pairsPThis 
is followed by appendixes which provide detailed deriva- 
tions of some of the technical results discussed in the 
main text. 



II. O(N) FIELD THEORY 

The relativistic O(N) field theory describes long- 
wavelength correlations of condensed matter systems 
with N real parameter fields with global O(N) sym- 
metry and two time derivatives in the LagrangianP For 
example, 0(3) theory describes unfrustrated Heisenberg 
antiferromagnetsp^ and 0(2) theory describes strongly 
interacting lattice bosons at commensurate^ and half- 
commensurate 5 fillings. The role of the speed of light is 
played by a spin wave or sound velocity. The Euclidean 




FIG. 2. Fluctuations in the Mexican hat potential. (€>) 

is the order parameter. The massive Higgs and massless Gold- 
stone modes can be represented either in terms of longitudinal 
(a) and transverse (ir) degrees of freedom, respectively, or in 
terms of scalar (p) and direction (<9 M n) degrees of freedom. 



time action reads as 

S = — ( d d+1 x 
2 5 J A 



(5,*) 2 + ^(|*| 2 -7V) 



(2) 



Here <& is dimensionless, the momenta integrals are cut 
off by the ultraviolet wave vector A, and the bare mass 
m is given by the microscopic scale of the system, of 
order A. 

Note that Eq. ([2| can be regarded as a coarse-grained 
version of the O(N) non linear a model. Very deep in- 
side the ordered phase, the a model does not allow for 
amplitude fluctuations; that is, the amplitude excitations 
are infinitely gapped. More generally, the model can be 
coarse-grained to a soft-spin model, as treated here, and 
the amplitude excitations have a finite energy. 

For d>2, there is an ordered SSB phase at weak cou- 
pling g < g c (A). Deep in the ordered phase (g <§; g c ) the 
order parameter is | | ~ y/~N. Quantum fluctuations 
reduce the order parameter to r(g, A)y/~N, while the am- 
plitude mode mass is reduced from its bare value m to a 
renormalized value m(g, A) until they both vanish at the 
quantum critical point g™. For Eq. ^ to describe the 
amplitude mode, the cutoff must be large enough to sat- 
isfy A ^> m, which holds particularly well in the vicinity 
of a quantum critical point. 

Fluctuations in the broken symmetry phase can be 
parametrized by 



(3) 



(rVN + a , tt) , 



where r is a constant, a is the longitudinal {not scalar!) 
fluctuation and 7r are N — 1 gapless transverse Gold- 
stone modes, as depicted in Fig. [2] Equation ^ can be 
expanded into harmonic, anharmonic, and counter term 
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parts, S — S + S A + S c 

S -1 

° ~ 2 3 JA 



d d+1 x 



(d^a) 2 +r 2 m 2 a 2 + (d^f 



S A 



S r 



d d+1 x 



25 

(r 2 — 1) m 2 

4g" 



r / , on 1 



(4) 



iV 



fT7T 



4A 

2r%/A <t + cr 2 + tt 2 



(cr 2 +7T 



2\2 



where the renormalized mass is m — rm , and r is de- 
termined by requiring (cr) = order by order in powers 
of g and/or 1/N. This counterterm prescription ensures 
that the Goldstone propagators remain massless at each 
order (see Appendix [A!) . 



x°= — n° = .; 
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FIG. 3. Infrared divergences at weak coupling. Dia- 
grams describing the infrared divergent weak coupling correc- 
tions to the longitudinal and scalar susceptibilities in Eqs. Q 
and ( II |. Solid lines are a propagators; dashed lines represent 



7r propagators. 



A. Dynamical susceptibilities 



Susceptibilities are generally defined by 



When the O(N) symmetry is broken, the generators 
T a fall into two classes: broken and unbroken. The A — I 



/i iaii mro rwo classes: DroKen ana unDroxen. j 
d d+\ x e n-x(A( x j B(0)) c , (5) broken generators rotate between $ x = a and $ 



where A — cr, tt, tt 2 , etc. Here the subscript c denotes 
the connected average, q — (ui n ,q) is the Euclidean mo- 
mentum, and ui n is a bosonic Matsubara frequency. The 
spectral function is XAsil^) = ' m Xab{1i w + *0 + ). 

In the broken symmetry phase, the longitudinal sus- 
ceptibility as defined by (I3J) is Xaa- ^ n contrast, scalar 
fluctuations are given by 



p{x) 



(\&(x)\ 2 - r 2 N) = 2ra 



N 



(6) 



One then defines the scalar susceptibility x ppJ which is 
related to a sum of cross-susceptibilities, 



\ op = r 2 (4x. 
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r 2 N X n 2 n 2 i 



(7) 



X n i„2 + 



r 2 N *~<j 2 a 



r 2 N X a 2^2- 



Anharmonic interactions S A in Eq. Q are responsible 
for cross-susceptibilities between the a and tt fields. In 
Sees. |IIl| and |IV| we calculate the longitudinal and scalar 
susceptibilities using weak coupling and large N expan- 
sions, respectively. Both approaches show that in the 
to — > (infrared) limit, xlaa an d Xs'ing are singular, while 
Xreg ls infrared regular. 

B. Conductivity 

For all N > 2, currents and conductivities can be de- 
fined as derivatives with respect to matrix gauge fields in- 
troduced into Eq. (B by setting <9 p * -> d^ + eA afl T a ^. 
The T a are the OyN) symmetry generators and e is the 
charge P 
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with j = 2, . . . , N. The remaining |(JV — 1)(A - 2) un- 
broken generators rotate among the (N — 1) components 
of the 7r field. In analogy with the 0(2) case, which has 
a single generator of the first kind (T — \cr v ), we define 
the (paramagnetic) currents of the broken generators by 



(8) 



their correlators as 

^(v') = (47WCM) - (9) 

and the generalized O(N) conductivity as 

cr( W ) = lm ^bb (w + j0 +, q = 0) - ^(|*| 2 ) 

5 (10) 

where the last term is the diamagnetic contribution. The 
O(N) conductivity is discussed more thoroughly in Ap- 
pendix Id] 



III. WEAK COUPLING LIMIT (g < 1) 

Diagrammatic perturbation theory for x m powers of 
g amounts to counting loops. To one loop order, 

XaM = X°(q) + X°(q)n (q)x°(q) + 0(g 3 ) (II) 

where x° = <?/(<? 2 + m ' 2 ) 1S the zeroth order longitudinal 
susceptibility and 77° is the polarization bubble, shown 
in Fig. [3j Since 77° is a convolution of two massless ir 



propagators, it diverges as (Appendix |B| 

jjO, , m jr>(N-l) f d d+1 k 1 

W 2N J (2n) D k 2 {k + q) 2 



N 



16\q\ 
32^2 



1 + ln 



(d = 2) 
(d = 3) 



(12) 



We ignore all other diagrams of the same order in g which 
do not contribute to the low-frequency dependence, such 
as the loop of massive er propagators. Analytically con- 
tinuing \q\ — > \J q 2 — (oj + ie) 2 and taking the imaginary 
part yields 



Xerer 



ng 



2v^ 



5(cj - *Jq 



(N-l) m y 0{uj 2 -q 2 ) 



32nN 

( 2tt 



(q 2 + m 2 



v 7 ^ 



1 



(d= 
(d= 



2) 
3) 



(13) 
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Note that Xo 

which is a direct consequence of the low-momentum di- 
vergence of II (q). This divergence is the quantum ver- 
sion of the divergent longitudinal susceptibility of the 
O(N) ferromagnet (N > 2) in d + 1 dimensions in its 
low-temperature ordered phased 

The leading order corrections of \ PP are given by 
the terms in Eq. Q which are dominated by infrared- 
divergent 77° factors. The dominant terms are depicted 
in Fig. [3| Computing the scalar susceptibility up to one 
loop yields 



Agr 2 



fq A r 2 



(q 2 + m 2 )' 



n°(q) + r 2 X (q) . (14) 



Note that the infrared singularities cancel out, leaving 
X pp to rise as q 4 II (q) at low q. As a consequence the 
scalar susceptibility decays rapidly at low frequencies, 




q 2 ) 3 / 2 Q(u ~ \q\) (d=2) 
q 2 ) 2 lnL-|<zl| (d = 3), 



(15) 



which enables a pronounced amplitude mode peak at 
uj ~ m due to the mass pole of %°. Note that 
when the action is parametrized in terms of p and n, 
the amplitude-direction coupling p(d^n) 2 gets two ex- 
tra derivatives relative to the longitudinal-transverse cou- 
pling (T7z 2 . This is responsib le fo r infr ared suppression by 
a factor lu 4 between Eqns. (13) and (15). This behavior 



is also seen at large N, as derived in Sec. IV and shown 
in Fig.[TJ 



A. Conductivity 

Diagrams contributing to a (to) are depicted in Fig. [4] 
The leading order conductivity (10) is of order g^ 1 and 



a) Oig- 1 ) b) 0(1) 




c) 0(5),'- 



d) 0(g) 





FIG. 4. Weak coupling expansion for the dynamical 
conductivity, (a) Order g -1 contribution to the weight of 
the superfluid delta function at to — 0. (b) Order g° diagram 
has a threshold at the amplitude mode mass. (c),(d) Two- 
loop self-energy and vertex corrections which contribute to 
the subgap conductivity. As N — » oo diagram (c) dominates 
and yields a ~ uj 2d ~ 3 . For N — 2 (the relativistic superfluid), 
cancellations between diagrams (c) and (d) suppress the sub- 



gap conductivity by four powers of lj, 
text). 



and 



(sec 



is all contained in the delta function weight at zero fre- 
quency, 

a{uj) = A6(lj) (16) 

where A = Ne 2 m d ~ 2 g~ 1 +O(g ). A nontrivial frequency 
dependence arises at order 0(1) from the ct-7t bubble 
diagram shown in Fig. [4j It exhibits a power law rise 
above a mass gap threshold: 



(7 (w) 



*S d e 2 
dw 2 m d ~ 



Airumi 



e( w 2 



(17) 



where S d is the surface area of a unit sphere in d dimen- 
sions. The 0(1) threshold conductivity is depicted in 
the main part of Fig. [5] Other 0(1) diagrams introduce 
frequency-dependent corrections at twice the mass gap, 
and overall renormalizations of the superfluid density A 
and the mass gap m. 

Two- loop diagrams are of order O(g). The two dia- 
grams which produce finite sub-gap conductivity are de- 
picted in the bottom of Fig. [4] As N — > oo, diagram (c) 
dominates the sub-gap conductivity. By power counting, 
the diagram scales as a ~ uj 2d ~ 3 , implying a significant 
sub-gap absorption for d — 2 for the large N conductiv- 
ity. 

At smaller TV, however, diagram (d) becomes compa- 
rable to (c) but opposite in sign, tending to cancel the 
sub-gap conductivity. After a lengthy calculation, shown 
in Appendix lEj we obtain the power series in frequency, 



~d=2 _ 

9 T \ 

zd— 3 ge^m 

9 ~ 37r 2 2 9 JV 



15m + 105rr; 

^ _ 
4m 3 T 10m 5 



32ui 



* +( 3Ar - 5 )3Tf^+- 



■}■ 



15m 5 

+( 9Ar - 16 )i^+-} • 

(18) 



Remarkably, we find that for N — 2, the coefficients of 
the two lowest powers vanish! This result can be under- 
stood as a consequence of the complete 0(2) symmetry 
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- g io- : 




layered, charged 
lattice bosons 



FIG. 5. Dynamical conductivity for 0(2) (relativistic 
bosons) in two dimensions. The arrow at zero frequency 
denotes is the superfiuid delta function response. For neutral 
bosons, there is a weak 0{g) sub-gap tail (see inset). The 
conductivity of a bosonic layered superconductor (red online) 
is plotted for using intralayer plasma frequency oj p = 10 m and 
a much smaller interlayer plasma frequency of ui^ = 0.1m. 



of the conductivity, which implies that it does not excite 
Goldstone fluctuations, similar to the scalar susceptibil- 
ity. In contrast for N > 2, the conductivity is not a pure 
scalar response, since it depends explicitly on the broken 
symmetry direction. 



IV. LARGE TV LIMIT 

The inverse number of components \/N controls an ex- 
pansion about the N — oo limit for all values of g. This 
allows us to approximate the finite N system in both or- 
dered and disordered phases, except close to the quantum 
phase transition at g c . Note that all anharmonic terms 
in Sfr are suppressed by negative powers of N. 

The large N renormalization of the order parameter is 
(see Appendix [A| 



r 2 (g,A) = l-g[ 



d d+1 q 1 =!_J_ 



(2n) d + 1 q 2 



(19) 



with r 2 vanishing linearly with g at the quantum critical 
point, 



4tt J A d = 2 
8n 2 /A 2 d = 3 



(20) 



The renormalized mass m 2 = w,q (l—g/g%°) also vanishes 
at 5 c°°- 

Now we evaluate the large N longitudinal susceptibil- 
ity, given by 



xZ(q) 



g 



q 2 + m 2 - gS a (q) 



(21) 



where E a is the longitudinal self-energy given by the RPA 
sunP, 



n°{ q ) 



1 + gn (q)/m 2 ' 



(22) 



as shown in Appendix |c| Note that, since 27° (q) diverges 



as q — > 0, S a (0) = rrr/g, in agreement with an exact 
Ward identitypES As a consequence, the pole in \ Ga at 
q 2 = m 2 gets replaced by a branch cut starting at q 2 = 0. 
For instance, for d = 2 we obtain 
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x<ya{q) 



16 |q|m 2 
\q\ 



(23) 



This infrared singularity, given at large N to all orders in 
g, agrees with the singularity obtained for all N at order 
g 2 in Eq. Q. 

Now we evaluate the scalar susceptibility [Eq. ([7])], 



Xpp 



^r 2 Xo 



sing 



0(1/N). 



and X T 2 T 2 are given by 



Ag 2 E <J {q)/m 2 
q 2 + m 2 - gScr (q) 



and 



r 2]\[ X-7r 2 7r 



'(«) = 



V 

-.2 _ 



+m 2 — gS a {q) 
Summing all the contributions in Eq. ^ yields 

4g r2 f, , sq 4 ^<7{q)/m i 



x P p{q) 



i 



(24) 



(25) 



(26) 



(27) 



+ \ q 2 + m 2 — gS (T {q) _ 

Note that the factor of q in the numerator suppresses the 
low q singularity in the denominator, and Xpp(q) °t \q\ d+1 
at low momenta, just as in the weak coupling case. In 
Fig. [Tithe large TV approximations for x'o<r{u) and Xpp( w ) 
are plotted for d = 2 . We take g/g™{A) = 0.84, inside 
the ordered phase. The amplitude mode peak is clearly 
visible in x'p P while it is difficult to detect in x'aa- 

A. Width of the scalar peak 

Thus far we have compared the scalar and longitudi- 
nal susceptibilities and have shown that a peak can be 
discerned in Xpp( w ) even in cases where it is hidden by 
infrared divergences in x'aa ( w ) ( see Fig- ■ Deep inside 
the ordered phase, g <C g c , the peak in Xpp( w ) becomes 
very sharp in relation to its energy oj = m. However, as 
one approaches the disordered phase, g — > g c , the relative 
width grows until close enough to the transition the peak 
in Xpp( w ) becomes broader than m. Thus, close enough 
to the transition, it becomes impossible to identify the 
Higgs energy from XppM- 
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We can study the width of the peak systematically in 
the large N limit. The scalar susceptibility at q = can 
be written as 



4.9 



l Im [F*(u)} 



ml | w 2 + ( w a - m 2)F(u)\ 2 ' 



where 



F(u) 



(28) 



(29) 



For example, from Eq. ( 12 ) we find that for d = 2 and 



N = oo, F(lj) — —iLo/{2^j), where 



32 ' 



Thus, 



4fl 



27a; 3 



ttIq (uj 2 — to 2 ) 2 + 47 2 w 2 ' 



(30) 



(31) 



Eq. (31 1 is peaked at uj — to with width 
;he other hand, for 7 > to, the width 



For 7 < 

smaller than to. On t : 
becomes larger than to and the peak is shifted to energies 
larger than to. Hence, 7 = to marks the point beyond 
which the Higgs mass cannot be determined accurately 
from the peak in Xpoi 1 ^)- 

Note that 7 grows linearly with g, whereas to = m^r 
vanishes at g — g c , according to Eq. (19 1. This is de- 
picted in Fig. [6] where it is seen that close enough to g c 
the width of the peak 7 exceeds the renormalized Higgs 
mass to. The question of how close one can get to the 
transition before this happens is nonuniversal, since it 
depends on the ratio A/to . For A = 2m , this occurs at 
g/g c — 0.96, corresponding to to/too = 0.19. For larger 
values of A/mo, the softening of the Higgs mode can be 
tracked to very low energies before its mass can no longer 
be detected reliably. 

For d = 3 and N — 00, we obtain 



where 



4.9 



?7 2 7rw 4 sign(cj) 
to 2 + r\ 2 ui 2 ( 1 + log ^2 



32tt 2 ' 



(32) 



(33) 



In this case, in contrast to d = 2, a sharp peak can be 
observed arbitrarily close to the critical point g c , that is, 
for arbitrarily small to/toq. 



V. EXPERIMENTAL PROBES 
A. Lattice bosons near the Mott transition 

Ultracold bosons in an optical lattice undergo a Mott 
transition at integer fillings as the lattice potential is 



g/gc = 0.1 



g/gc = 0.5 



X 



g/gc = 0.9 
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FIG. 6. Dependence of the scalar peak on g. Scalar sus- 
ceptibility x"p( w ) f° r different values of g/g c in d — 2 dimen- 
sions. Results shown in the large N limit for A = 2mo. (In- 
set) Renormalized mass m (solid curve) and width 7 (dashed 
curve) of the peak in Xpp( w ) as a function of the tuning pa- 
rameter g, expressed in units of the bare mass mo. For g — > 0, 
deep inside the ordered phase, the mass approaches its bare 
value mo and the peak is very sharp, 7 — > 0. As g ap- 
proaches g c , the mass is softened and the width grows, until 
for g/g c > 0.96 the peak energy is smaller than the width. 



strengthened.^ The transition is well described by an 
0(2) relativistic field theory,^ with the lattice strength 
controlling the radius of the Mexican hat. Amplitude 
oscillations could, in principle, be observed after quench- 
ing the Mott phase into the superfhiid phase, similar to 
coherence peaks recovery seen by Greiner et al.p^ who 
studied the opposite quench direction. More directly, 
by modulating the lattice potential at some frequency 
uj, one can excite the scalar mode of the superfluidP^ 
Thus, the system is predicted to absorb energy at a rate 
wXppH- Such measurements have indeed been carried 
outp^l and although no sharp peaks were observed, this 
can likely be attributed to the presence of a harmonic 
trapping potential, whose effect is to smear the energy of 
the amplitude mode, and the application of long-lasting 
modulation pulses that pushed the system outside of the 
linear response regime. 

The phase of the optical lattice can also be modulated 
instead of its amplitude. The energy absorption rate is 
then given the optical conductivity cr(u>) at the phase 
modulation frequency^ which, as shown above, has a 
threshold at the Higgs mass. 

Finally, recent experiments using Bragg spectroscopy^ 
have shown evidence of the amplitude mode in interact- 
ing lattice bosons. Although these experiments were car- 
ried out beyond the linear response regime treated here, 
they demonstrate that the Higgs mode is, in principle, 
observable though the use of these probes. 
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B. Raman scattering in Antiferromagnets and 
Charge Density Wave Systems 

Heisenberg antiferromagnets, and incommensurate 
charge density wave systems can be effectively described 
by relativistic 0(3) a nd Q( 2) theories, respectively. In- 
elastic light scattering * 20 * 21 ! can effectively couple to the 
square of the local order parameter. For example, in 
antiferro magne tic insulators, light couples to bond spin 
operators, 1 22 1 23 1 



^Raman ~ D / ! E tl E 'n S X ' S x+r, 

X,T) 

cx^T jd d xE v E' n (\*{x)f 



(34) 



Thus, with incoming (outgoing) electric field polarization 
E (E 1 ), consistent with the crystal symmetry, the Raman 
spectrum measures the scalar susceptibility Xp P (w). The 
small but finite separation between the two spin oper- 
ators leads to small corrections, denoted by ... in Eq. 
(34 1. The most important of these corrections, propor- 
tional to (17 • V7r) 2 gives a direct coupling of light to 
pairs of spin waves and leads to a broad background sig- 
nal. The background is well-behaved in the infrared due 
to the extra spatial derivatives in the coupling and hence 
is not expected to hide any peaks present in XpoO^)- 

Indeed, a pronounced Raman peak has been measured 
in several magnetic compoundsp2 Theoretically, it has 
been analyzed as a two-magnon resonance, or bound 
statej^ which is equivalent (in its quantum numbers) 
to the amplitude mode. We note that inelastic neutron 
scattering probes the longitudinal susceptibility of anti- 
ferromagnets, since the neutron spin couples locally to 
the Neel vector. At the Bragg wave vectors, the ampli- 
tude mode peak is therefore expected to be obscured by 
the singular low-energy scattering. 

For incommensurate charge density wave the order 
parameter has a massive amplitude mode, and a soft 
translational mode. Light excites the amplitude mode 
by inducing inter-band transitions This has been used 
to detect the amplitude mode in low-dimensional CDW 
systems through Raman scattering^ an( j femtosecond 
pump-probe spectroscop y 2 " 29 * experiments. In contrast 
to antiferromagnets, in the case of CDW order, neutrons 
couple to the local charge density and hence act as a 
scalar measurement. Indeed, the amplitude mode of a 
CDW has been measured using neutrons Similarly, 
neutrons have been proposed to detect the amplitude 
mode of a DDW stateP 



C. Superconductors 

Granular superconducting films and low-capacitance 
Josephson junction arrays which exhibit superconductor 
to insulator transitions, can be effectively described in 
terms of a bosonic 0(2) relativistic field theoryP3* By 



proximity to a Mott insulating phase, the amplitude 
mode may be suppressed below the BCS pairing gap, 
and appear as a long-lived collective excitation. 

In homogenous BCS superconductors, the detection 
of the amplitude mode by Raman scattering has been 
proposed in the presence of a coexistent charge density 
waveP^l Here we propose that the amplitude mode may 
be observed in a class of "bosonic" superconductors, that 
is, those exhibiting short coherence length, low super- 
fluid de nsity, and perhaps a pairing gap above T c (e.g., in 
cuprates^El). Such superconductors may be described 
by charged lattice bosons, which may be treated by the 
0(2) theory of Eq. ^ with Ion g-rangc Coulomb interac- 
tions. The optical conductivity of HCB couples directly 
to the amplitude modej^ however, since Cooper pairs 
are charged, long ranged Coulomb interactions modify 
the Goldstone mode's dispersion. In a three-dimensional 
sample, the phase fluctuations are gapped at the plasma 
frequency uj p . 

Nevertheless, in highly anisotropic layered supercon- 
ductors such as the cupratesp^ the threshold for opti- 
cal absorption at zero temperature is shifted from m by 
the relatively small c-axis plasma frequency w|j«Up, as 
shown in Appendix [FJ We propose that Raman scatter- 
ing, which couples to the 0(2) scalar susceptibility, may 
be used to observe the amplitude mode peak. The mass 
m of this mode is expected to be of the order of the su- 
perfluid density and T c and decrease toward the quantum 
phase transition into the insulating phase. 

For example, in Bi 2 Sr 2 CaCu208+5, w a & ~ leV, while 
lo c ~ I meV, leading to conductivity qualitatively similar 
to the charged latice bosons shown in Fig. [5} We propose 
that the amplitude mode may be partially responsible for 
the rise in optical conductivity in the mid infrared regime, 
as observed above 400 cm -1 El. 



VI. SUMMARY 

We have calculated the scalar susceptibility and the 
conductivity by weak coupling large N expansions within 
the SSB phase of the two- and three-dimensional O(N) 
field theory. In contrast to the longitudinal susceptibil- 
ity, where low-frequency dissipation by Goldstone modes 
can entirely mask the amplitude mode, the O(N) sym- 
metric susceptibility exposes the amplitude mode as a 
finite- width peak that is uncontaminated by infrared con- 
tributions arising from the Goldstone modes. Similarly, 
the 0(2) conductivity exposes the amplitude mode as a 
broadened threshold in frequency. 

From an operational point of view, the mass m can 
be extracted experimentally from the peak in the O(N) 
susceptibility. While the determination of the peak en- 
ergy is limited by the width of the peak, one can identify 
such a peak unambiguously as arising from a Higgs mode 
by tracking the peak position m as a function of a tun- 
ing parameter near a quantum phase transition. Deep 
in the ordered phase, the peak is very sharp relative to 
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its energy m. As one approaches the phase transition, 
the peak energy softens and its width grows. Tracking 
the softening of the peak position is then possible except 
for a region very close to the transition, when the peak 
width becomes comparable to its energy. Similar consid- 
erations apply to the broadened threshold in the 0(2) 
conductivity. 

The suppression of sub-gap absorption by four powers 
of frequency was derived by precise cancellations between 
self-energy and vertex corrections in the conductivity. It 
is easier to understand in the amplitude-direction rep- 
resentation, where p is coupled to the derivatives of the 
order parameter direction p(d^n) 2 . Thus, the intrinsic 
dissipation of the p self-energy is suppressed by four pow- 
ers of momenta. The current operators in all non-Abelian 
theories N > 2 are not O(N) symmetric, since they con- 
tain explicit coupling to the angle variables. For 0(2), 
however, the current operator j x = epd x ip is rotationally 
invariant and therefore has four higher powers of sub-gap 
absorption than the non-Abelian conductivities. 

Our conclusion is therefore that the amplitude mode 
is in fact long-lived even for moderate g < g c , but its 
detection requires using 0(N) symmetric experimental 
probes. We propose Raman scattering for antiferromag- 
nets and superconductors; coherence peak oscillations 
and lattice modulation experiments in supcrfluids near 
the Mott transition; and optical conductivity in bosonic 
superconductors. 

Acknowledgements. We thank Ehud Altman, Nicolas 
Dupuis, Rudi Hackl, Sebastian Huber, Steve Kivelson, 
Ancesh Manohar, Subir Sachdcv, and Wilhclm Zwerger 
for useful discussions. Wc acknowledge support from the 
Israeli Science Foundation, the U.S. -Israel Binational Sci- 
ence Foundation, and the European Union grant IRG- 
276923. One of the authors (DPA) was also supported 
by NSF grant DMR-1007028. We are grateful to KITP 
at Santa Barbara and Aspen Center For Physics where 
some of this work was initiated. 



Appendix A: Counterterms in the broken symmetry 

phase 

At small g < g c , the order parameter acquires a finite 
vacuum expectation value (VEV) (<&<t) = r\/~N in the 
a direction. We choose to expand the Euclidean-time 
action about the VEV 



* = (rVN + <t,tv) 
which leads to the action 
S = S Q + S A + S c 



(Al) 



1 

25 



d+l 



(d^af +r 2 ml<7 2 + (£> M 7T) 



(A2) 



25 J a 



d d+l 



(- 3 + - 2 ) + T^(- 2 +- 2 ) 



g c= (7 % 1)m ° l<r' +i - 



AN 

2rVN a + a 2 + tt 2 



The harmonic action 5 provides zeroth-order massive 
and massless propagators, 



G L = 



g° = 



k 2 + in 2 

5_ r 

k 2 0ij ' 



(A3) 



where m — r mo is the renormalizcd mass. 

The parameter r is chosen such that that the VEV of a 
is zero, that is to say, that the fields a and 7r are expanded 
about one of the true ground states of the system. This 
is equivalent to requiring the vanishing of the a tadpole; 
that is, the sum of all 1PI diagrams with a single external 
a line must vanish. 

In the large TV limit, we can compute r in closed form. 
At leading order in N, 0(y/N), there are only two 1PI 
diagrams: a sigma line terminating in a 7r loop, 



mgr(JV-l) 



d d+1 k g 
2gVN J (2^) d+1 fc 2 "' 

and a a line ending in a counter term vertex, 

m%r(r 2 - l)y/N 
2^ • 



(A4) 



(A5) 



Setting the sum of the two terms to zero yields (for N 

oo), 



r 2 = 1 - 5 



/ 



d d+1 k 1 



+ o 



(2iT) d + 1 k 2 \N / 

= l-5/5c°°- (A6) 

Evaluating the integral with a cutoff A on the spatial 
momenta (but no cutoff on the Matsubara frequencies) 
yields 



[a-k/A d = 2, 
9c \8n 2 /A 2 d = 3. 



(A7) 



Since we expand the fluctuations about the true ground 
state, Goldstone's theorem guarantees that 7r field is 
massless. Indeed, in the large limit, the 7r self-energy 
is the sum of two diagrams, a 7r loop and a counterterm: 



= -6-^ 
13 25 

= 0. 



d d+1 k g 

(2TT) d + 1 k 2 



+ (r 2 1) 



+ o 



1 

TV 
(A8) 



Thus, for A^ = oo, the bare and renormalized 7r propa- 
gators are identical to each other. 

A similar computation shows that the constant 
(momentum-independent) contribution to the a self- 
energy cancels for N — oo, thus identifying m = m r 
as the renormalized mass of the amplitude mode beyond 
tree level. 

In the weak coupling regime, all of our computations 
are carried out to leading nontrivial order in g. In this 
case, we can set r = 1 and also ignore the counterterms, 
which only correct our results at subleading orders in g. 
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Appendix B: Computing the longitudinal 
susceptibility in the weak coupling limit g <C 1 



In momentum space, the bare susceptibility Xa-aio) is 
XIM = » (Bl) 



where q 2 = q^q^- The full susceptibility is 
Expanding, we have 



(B2) 



x„M = Ah-2 + \Ah^) K(q) + o(g 3 ) ■ (B3) 

q z + m z \q + to / 

For small q, is dominated by the polarization in- 

sertion n°(g), computed in the next section. We write 

E (T (q)=n°(q) + ... (B4) 

where . . . denotes terms that are either of higher order 
in g or infrared finite. 



1. Polarization insertion 



The integral we must do is 



N 



(<?), (B5) 



Here D — d + 1, where d is the spatial dimension, and 

= 1 f d D k 1 
D(q> ~ 2 J (2n) D k 2 (k + q) 2 

A/q 



(B6) 



q D ~ A f, n_o /" 



2(2tt 



, W 3 ' 



o 



p 2 + 2p cos 9 1 + l 



The D-dimensional unit vector is 



n = I cos 6^ , sin 6^ cos 6> 2 , ■ • ■ , sin • • ■ sin u D _ 2 cos <p) , 

and the metric is 

d(2 D = (sm D - 2 6 1 d6 x )--- (sin 9 D _ 2 dO D _ 2 ) d<j) (B7) 

Here D = d + 1, where <i is the spatial dimension. 
For <i = 2, we have D = 3 and we can take A — > oo, 

OO 71 271 

h(q)^ 2 L q ^fdpfd9sm9fd ( f ) 



p 2 +2pcos8+l 





oo 1 



p 2 +2 i a; p+i 



-1 

oo 1 

1 / dp J n P+l _ 1 / d£ m 1+p 
&7T 2 q I p p— 1 4tt 2 9 / p 1— P 


2^(1+3^ + ' 



16 q 



For d = 3, we have D — 4. We must retain the ul- 
traviolet cutoff A, which we take for convenience to be 
isotropic in the spatial and temporal dimensions. We 
make use of 



dtp 



1 



2tt 



We have 



a + b cos ip ^Ja 2 - b 2 



A/q 



Q(a 2 - b 2 ) . (B8) 



7 4(<?) = 1(2^)4 Jdpp Jd9 sm 2 9 ■ Att ■ p2+2p l osg+1 
o o 

A/q tv 
1^3 p ?+2pcosfl+ I 

-7T 

4/g 



= s^ dp 





A/9 



3^2 /dp 



(p 2 +l)-|p 2 -l| 
4p 



(p 2 +l) 2 


i 






|p 2 -i|_ 




i 




//1 2 \1 




1 4- In ^ 


32tt 2 




\r )\ 



(B9) 



These results for d = 2 and d = 3 are consistent with 
the general results in Eqs. (13) and (14) of Ref.^. 



2. Analytic continuation to real frequency 

We set 

q = \Jq 2 - (u) + ie) 2 



\/q 2 - u 2 



if lu 2 < q 2 
/ 2 \/uj 2 - q 2 if uj 2 > q 2 



(BIO) 



The susceptibility is 



Xaa q 2 + l-{uj + ie) 2 

iJ2 



, (N - l)mjr 



1 



16 AT [q 2 + 1 - {lu + ie) 2 } 
(d = 2) 



1 

2^ 



1 + ln 



(q 2 -(w+ie) 2 ) (d— 3) 



(Bll) 



Therefore, 
a _ 



ng 



5 (to — \J q 2 + ra 2 ) 



(N-l) m y e{co 2 -q 2 ) 



32nN {q 2 + m 2 - lo 2 ) 2 
I 1 (d = 3) 



(B12) 
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+ 



+ 



w iy ^H^r- + 



+ 



The RPA captures the leading order behavior in the N — > 
oo limit. Thus, 



q 2 + m 2 - g^ PA (q)' 



(C2) 



Similarly, the scalar susceptibility is obtained from 
Eq. using the RPA sums shown in Fig. [7j 



v = «'E * + •£ * — • £ * + 



Xr°°(9) 



1H g g 4 sg FA (g)M 4 

g 2 + 77i 2 I q 2 + m 2 — gY^ FA {q) 



4gr 2 



(C3) 



FIG. 7. Diagrams in the large- N limit. For N = oo, only the 
RPA sums shown here contribute. 



Appendix D: Conductivity of O(N) models 



Appendix C: RPA and large N theory 

Let us evaluate 'E (T (q) to order K in perturbation the- 
ory, where we include two air 2 vertices and K (ir 2 ) 2 
vertices (see Fig. [7]). There is an overall factor of 
(— l) K+2 /(K + 2)1 from the exponential. Selecting two 
of the cr7T 2 vertices from the K + 2 terms results in a 
combinatoric factor ( if ^ 2 ) ■ The two air 2 vertices can be 
interchanged, and the K (7T 2 ) 2 vertices can be permuted, 
resulting in a factor of 2! • Kl. When we contract one of 
the 7T legs of the first air 2 vertex with the first (w 2 ) 2 ver- 
tex, there are four choices of legs from the latter vertex to 
choose from. For the second leg, there are three choices, 
but if we want to maximize powers of N there is only 
one choice. (Recall that we are contracting a-K a -K a with 
7r^7r /3 7r 7 7r 7 .) We are left with two uncontr acted legs of 
the first (tt 2 ) 2 vertex, and we get another factor of four 
from the second (w 2 ) 2 vertex. After running through all 
K of these (7r 2 ) 2 vertices, we have two remaining legs 
to contract with the 7r legs from the second an 2 vertex, 
yielding two possibilities. Each of the K (tt 2 ) 2 vertices 
comes with a factor rriQ/8gN, and each of the two <ttv 2 
vertices comes with a factor rri^r /2gy/~N ' . After all the 
legs are contracted, we are left with [K + 1) loops, each 
containing a n a propagator at momentum k + q and a 
7r Q propagator at momentum —k. The propagator at 
momentum k is g/k 2 , where k 2 = k^k^. For each of the 
loops, there are N — 1 choices of the vector index a for 
TT a . Putting this all together, we obtain a contribution 



l^lll . ( K - 1 V 2! • K\ ■ A K ■ 2 • [0^ ■ f^f 



(K+2)[ I 2 



\8gN J 



(/ 



A'+l 



K 



Summing this over all non-negative K, we obtain the 
RPA self-energy for the a field, 



n°( q ) 



l+gn°(q)/m 2 



(CI) 



An O(N) field theory can couple to a set of gauge, 
fields, with which currents and conductivities can be de- 
fined. Consider a local gauge transformation 

-> 0{x) #(ar) = e 9 ^ x) T " , (Dl) 

where {T a } are the ^N(N — 1) generators of the group 
O(N). The generators are real antisymmetric matri- 
ces which we normalize according to the convention 
j r (T a T b ) = -2S ab . A convenient basis then is the set of 
matrices 



(D2) 



where a denotes the composite index (/, J), where I < J, 
which runs from 1 to ^N(N — 1). For this basis we have 



(D3) 



The gauged O(N) model is defined by the Lagrangian 
density 

£e = i (0 M * + A^f + ^ (|*| 2 - Nf , (D4) 

where is an antisymmetric tensor vector potential 
which can be expanded in the generators, viz., A — 
A a a T a ■ Gauge invariance follows from the gauge trans- 
formation rules, 



A„ -> T A„0 - T d„ 



(D5) 
(D6) 



There are ^N(N — 1) O(N) currents, one for each gen- 
erator. We have 



SA al j,(x) 
1 

I" 



1 



9 ' 9 



(D7) 



(D8) 
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where Iujx is the paramagnetic (diamagnetic) current. 
The corresponding Kubo formula is 

(/ aM (x)> = - jd D x'K^(x,x')A bv {x') + 0{A 2 ) , (D9) 



with 



Kf v {x,x') = (l ail {x)I bv (x 1 )) 



(D10) 
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5^5\x-x')(T a <S>(x)-T h 3>{x)) 



We separate if into paramagnetic and diamagnetic con- 
tributions, with Kf v = K^ b + and K^ b (x, x') = 

<C(z)^(z')>- 



1. Symmetric phase 

In the symmetric phase, (3?) = 0, and the response 
function K^ b (x, x') is diagonal in (a, b). Summing over 
all O(N) generator indices, we define 



JjV(iV-l) 



Ka V {x,x') 



N(N ■ 



— K™{x,x'), 



which equals 
1 

N(N - l)g 2 
2 

~Ng 



V*i u ^] W 3 U ^i)x X*i U i>^3 *J L, "*Vx' 



S^(x-x')(^( X )^(x')). 



(Dll) 



Note that we have normalized by dividing by the total 
number of generators. We can define an O(N) conduc- 
tivity a {to) as 



cr(w) 



uid 



(D12) 



/x=l 



where the sum is over the spatial values of the space-time 
indices. 

Writing $ = (N + p) 1 ^ 2 n as before, we have 



d^ = {N + p) 1 ' 2 d il n- 



2(N + p) 1 / 2 



Making use of Eq. (D3), we can write 



(D13) 



(D14) 



2 
Ng 



S^ix-x') ((N + p) x (N + p) xl ) 



where 



Ji 



■A x ) = {N + p) x (n i d tl n J - n^n^ 



(D15) 



Note that for N — 2 we have n = (cos if , sin cp) and 
(n 2 d li n 1 — n 1 d^n 2 ) — d^ip, in which case the above 
expression reduces to a familiar form. 



2. Broken symmetry phase 

When the O(N) symmetry is broken, the generators T a 
fall into two classes. We define class A generators as those 
which rotate between <S> X = a and <P 1+: > = Wy There are 
(N — 1) generators of this class, with j = 1, ...,N — 
1. Class B generators rotate between <P 1+: > and . 
There are | (N — 1)(N — 2) generators of this class. Note 
that the total number of generators in classes A and B is 
(iV-l) + i(iV-l)(iV-2) = \N{N — 1), the dimension 
of 0{N). Thus, we can take 



T u = S kA,i ~ 5 iJk,i (Ki) 



(D16) 



T, 



S kAf ~ 6 U S kJ> (K i < i') ■ (Di7) 



The response function K°" v is diagonal in the generator 
indices, so we can in principle study two response func- 
tions, Kf£ and K**. Note that = K*£ = 0. 
For the class A generator T a with a = (j), we have 



IL = 1 (^d ll a-(rVN + a) d,,ir. 



a M g \ "J "f 



(D18) 



The diamagnetic contribution to the response function is 
K^\x,x') = --J ab 5^(x-x') 

'(rVN + a(x)) 2 +7rJ(a;)) . (D19) 
For the class B generator T a with a = (jj 1 ), we have 



(D20) 



The diamagnetic contribution to the response function is 

K^ b (x,x>) = -h^J^ix - x')(it%x) + 4(x)). 

(D21) 

Note that for N = 2, class B is the empty set. Hence, 
in analogy with the 0(2) conductivity, in what follows 
we focus on the class A response function and will drop 
the AA superscript. Then, averaging over the (N — 1) 
generators in this class, we have 



(D22) 



K^(x,x') = {N r l 1)g 2 (d^TT x ■ d„TT 



r z^jg2 (d^ x ■ (ad u ir - Trd u a) x , 



{N-l)g* 



(Jv-i)s 2 



[ad tt - tt8 a) ■ {od v ^ - ird u a) 



We are interested in the imaginary conductivity at 
q = and finite frequency, which we compute up to two- 
loop order. Thus we can omit the diamagnetic term, and 
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the first three terms of Eq. (D22). Taking the Fourier 
transform, we then obtain 



1 



d D x e iq-(x-x') 



(N - l)g 2 



(D23) 



where q ■ x = q^x^ . In what follows, we valuate this 
expression to one- and two-loop order and use Eq. (D12) 
to compute the optical conductivity. 



Appendix E: Optical conductivity to order g 
1. Conductivity at order g° 



At order g°, Eq. (D23l factorizes into the product of a 



a propagator and a n propagator. This yields the one- 
loop integral, 



d d k (2fc^ + q^)(2k v + q u ) 



x) d (k 2 + l)(k + q) 2 
We set fi = v = x and q = (iuJ n , q = 0). We then have 

K (i, , \ — 1 \ ~" f d d k fc 2 



(El) 



2 k 2 +{iu m +iuj n y 2 



1\ 

k A [dkk d+1 iV — 1 (E2) 



(A) *--4- 




FIG. 8. Order g diagrams contributing to the optical con- 
ductivity for uj < 2m. Solid brown lines are a propagators 
while dashed blue lines are tt propagators. Note that dia- 
gram (c) does not contribute to the subgap conductivity and 
is therefore not shown in Fig. [1] 



and the conductivity at order g° is 



<x(w) = Im K xx (u) + ie,q = 0) 



(E6) 



dk 



dui {2ir) d J """ ^W+l 




s(u-k-y/k 2 + l) 

-s(u> + k+\/k 2 + l\ 



where a k = \/k 2 + 1 and b k = k. 

Now the bosonic Matsubara sum can be written as 

i ]T J = - E »(P)Res[MP)] > (E3) 



where the sum is over the poles v of /i, and where n(y) = 
[exp(z//T) — l] is the Bose function. Note that n(±6— 
iw n ) = n(±b) when uj n is a bosonic Matsbara frequency. 
Hence, 



F(a,b,iuj n ) = 



I " ( ' a) r ^ r (E4) 

2a[(iu„-a) 2 -b 2 \ 2a[(iui n +a) 2 -b 2 \ 



+ 



n(-b) 



»(fe) 



21, 



[(iui n +b) 2 -a 2 ] 2b[{iio n -b) 2 -a 2 



Thus, at T = 0, 
F(a fc ,6 fe ,w + 'ie) 



a;-t-ie+a, +b. uj+ie—a.—b. 



(E5) 



2. Conductivity at order g 

There are 11 two- loop diagrams which enter the con- 
ductivity at order g. Of these, most either renormal- 
ize the zero-frequency superfluid stiffness peak or have 
a threshold at high frequency. This leaves three contri- 
butions, diagrammatically represented in Fig. |HJ which 
affect the finite frequency response below uj = 2m. We 
have 

q = (Q,iu j ), k = {k,iv m ), p = (p , ifj . (E7) 

We find it convenient to define 

fc 1 = fc, fc 2 = — P, fc 3 = P — fc, (E8) 

so that fc-L + fc 2 + fc 3 = 0, i.e. the vectors k l 2 3 form the 
legs of a triangle. 

The contributions A A , K B , and A',-, are given by 

K - —k i (\r 4 4 f ddkl f ddk2 

n A — (iV-l)g 2 8JVg 2 W v d I ( 27T ) d / ( 2w ) d 

P ^ 1 [k 2 +m 2 -(iv m+iu> .) 2 Y k 2 -{iv m ) 2 

. 3. 9 C E n\ 

fc!-(i?„+iw,-) 2 fc§-(i«„-^ m ) 2 ' V ; 
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K B — {N 



and 
K, 



1 1 ( Kf -\\ O 4 / g fc l / g fc 2 

{N~l)g 2 SNg 2 ^ JV L > ° d (27r) d I (2n) d 



^ k 2 +m 2 — (iv m +iu)j) 2 k^+m 2 — (i£ n ) 2 

* 9 - (E10) 



k 2 -(tv m ) 2 kl-{iv m -iU) 2 fc 2 -«„) 2 ' 



C — (N-l)g 2 SNg 2 



d d k 1 


f ddk 2 







a E g fc 1 ■ fc 2 



fc 2 -(ii/ m ) 2 /c 2 +m 2 -(iw J .+iC„) 2 
2 g 



fc 2 +m 2 -(w m -i£„) 2 fe 2 +m 2 -(w,„+it,; J .) 2 fe|-(i?„) 2 ' 

(Ell) 



Carrying out the Matsubara sum using Eq. ([E3j), first over £ n and then over z/ m , we obtain 



g(N-l) f d% f d% 



E Res 



fcl n(v) 1 1 

k 2 { V+ iu. + k 2 y-kl V 2 -k{ [ (j/+i ^)2_ fe 2_ m 2 



fc? n(v) 1 

k 3 (,u+iUj-k 3 )*-kl v2 ~ k 't [(^+^) 2 -fc 2 -m 2 ] 2 



R- - 9 [ d " kl ( d " k2 (k 2 I A- 2 fr 2N l V Rr-. 

IX B — Nd (27r) d / (27r) d ^1 + K 2 K 3J KeS 



(E12) 



V fe 2 2 +™ 2 ( l ^- v /fei+™ 2 ) 2 -fc| (<^+ N /fcI+m 2 )-fel ^- fc ? (^+^) 2 -fc 2 -r 



and 



, n(t/) 1 1 1 1 

k 2 {iu>.+k 2 ) 2 -k%-m 2 (v+iuij + k^-k 2 v 2 -k 2 (v +iul .) 2 -k 2 - 

, n(v) 1 1 1 1 

fc 3 (;y-fe 3 ) 2 -fe 2 (iy+iw 3 -fe 3 ) 2 -fc 2 i^-fc 2 (i/+iw j ) 2 -fc 2 -m 2 

K _ 3g f d% f d% ,, 2 ,2 ; 2^n J n(^) 
- Nd (27r)rf / (2tt)<< ^1 + K 2 K 3 J 2^ | ^ 2 -fc 2 



(i^+ic^) 2 — fc 2 — m 2 
1 



fc 2 (i/+k 2 ) — fc 2 — m 2 fiuj- — fe 2 l — A; 2 ,— m 2 



V / fc|+™ 2 (W^ 3 +^fc 2 +m 2 ) -fc 2 -m 2 (ii^+^fcf+m 2 ) -fc 2 

1 1 1 

v/fcf-H™ 2 (i^-^/fcl+m 2 ) 2 -*: 2 ( y - k 2 +m 2 ^ - k 2 - rr. 



Using Eq. (D12|, we next obtain the 0{g) conductivity from these expressions. 



a. The O(g) contribution ita(w) 



In computing the contribution cr A (w), we must compute 



E Res 



= E Res 



1 

46? 



m 

(u + iu} j + by) 2 [y + iwj - by) 2 
F{-iuo j - by) + ^ F'{-iui j - b{) - F{-iu j + by) + by F'i-iu), + by) 2 



(E13) 



E 



- [{ V + ^u: 1 ) 2 -b 2 y]■ 



Res[F(i 
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where b- = ^Jk 2 + m 2 , with j <E {1,2,3}. Note that at T = 0, n(—iujj — b x ) = —1 and n{—iuj- + 6 1 ) = 0, where 
n(z) = (cxp(z) — 1) _1 is the Bose function. We now set iuj n = u + i0 + , in which case 



„ (ui\ = g^- 1 ) im i /" d fcl /" d fc 2 

AV ; W J (2n) d J (27r) d I 46ffe 2 fe 3 ^(w+b 1 ) 2 -fcf 1 [(^t^-fc?] 2 ; ^ft-t 2 ) 2 -^ 1 (b 1+ k 3 r-ki 



J |_ 26 1 (gj+t» 1 ) 



_&3 | &2_ 



, _Jn 1 I 2(b 1 -fc 2 )fc 3 2(b 1 +fc 3 )fc 2 

4b?fe 2 fc 3 (^+fc 1 ) 2 -fe? U( b i- fc 2) 2 -fe|] 2 [(6i+* 3 ) a -*2] : 



fc 2 



Combining terms, we find 



g(^-i) 

Nd 





2fe 2 


d% 


/"d d fc 2 




(2^) d 



2fc 2 fe 3 [ (w _ fel )2_ 6 2]2 ^ (w _ /Si+fe2 )2_ fe 2 (u ,_ fei _ fe3) 2_ fe 2 

kl 1 1 \ 

2fe 2 fc 3 [(fe 2+ fe 3 )2_ 6 2]2 ( aJ + fc2+ fc 3 )2_fc2 C ' 



J |_ 2t» 1 (fa;+i> 1 ) 



(fe 2 +fe 3 )^( w +6 1 ) 2 -fe? [(uM-b^-fc?]' 



(E14) 



26 1 fe 2 fe 3 ( W +fe 1 )2_ fc 2 [ fc 2_ (fe2 + fe3) 2] 2 2 fe 2 fe 3 [(cu-fc^-fc 2 ] 2 (w _ fei) 2_ (fe2 + fe3) 2 

fc? 1 1 I 



2/c 



,fe 3 [(fc 2 +fc 3 ) 2 -6 2 ] 2 (cu+fc 2 + fc 3 ) 2 -fc 2 ( ■ 



(E15) 



Now we use 



Im 



[(w - /c) 2 - fe 2 ] 2 46 2 



[J'(w - 6 - fc) + 5'{u + b - k)] 



(E16) 



and set oj > to obtain 



a A {u > 0) 



ng(N-l) r d% r d% 

J (27Tpy (27Tp 



2fcj <5(u; — fc 1 — fc 2 — fc 3 ) &i(& 2 + ^3) ~ ^i ~ &i) 



fc 2 fc 3 



b\ - (k 2 + fc 3 ) 2 



6 2 fc 2 fc 3 6 2 (fc 2 + fc 3 ) 2 



(E17) 



The first term inside the large braces integrates to a result which is finite for arbitrarily small oj. The second term 
yields a threshold behavior and is finite only for oj > m. 



b. The O(g) contribution ctb 

Summing the residues, we find 



JLUl / _ fc l / _ "2 ^1.2 _i_ 1,2 _ j,2~\ J _1 1 1 1 

k^b 2 (w-b 2 ) 2 -fc 2 (6 2 -fe 1 ) 2 -fe| (w-fc 1 ) 2 -6 2 

+ -J I I 

fc i fe 3 (fci+fe 3 ) 2 - fc l ("-fci-^a) 2 -*:! (' 
.1 1 1 1 1 1 



1 


1 


1 




1 


k 1 k 2 
1 


(cu+fc 2 ) 2 -6 2 
1 


(u-k 1+ k 2 y 

1 




(^-fe 1 ) 2 -6 2 
1 


b l b 2 
1 


(w-fe 2 ) 2 -fc| 
1 


(cu+b^b,) 2 - 

1 


it 2 


(uj+b^-kf 
1 


b l fe 3 


(w+bj+fcjj) 2 - 


-fc 2 (b 1+ k 3 r 


—k 2 


(w+b^-fc 2 


1 


1 


1 




1 


k 2 k 3 


(a) + fe 2 )2_ b 2 


(^+fe 2 +fc 3 ) 2 - 


-k\ 


(fe 2 +fe 3 ) 2 -b 2 



I I. I. 1 



(E18) 
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After some grinding, we obtain 

*b(" ><>) = - A 



rl d h I" rl d h 



b,k,k 2 (b 1+ fc 2 )2_ fe 2 (bi+fei+fe2 )2_ b 2 



(2ir) d J (2tt) 

1 1 1 I __L 1 1 

KKk 2 (6 1 -fc 2 )2-fc2 (6 1 +fe 1 -fe 2 )2- b 2 I" kibib2 (b2 _ ki) 2_ k 2 ( fe i+hi _ b2 )2_ fe 2 

1 1 1 i __1 1 1 

fc^b, (6 2 + fe 1 )2-fc 3 2 (b 1+ b 2 + k 2 y-kj + fc^fe, (fc 1+ fe 3 )2_ 6 | (bi _ fe3 )2_ fe 2 

8{uj — b 1 — fc-J 



k 1 b 1 k 3 (k ± — fc 3 ) 2 — b\ (fc 1 +fc 3 ) 2 — k 2 



b 1 b 2 fc 3 (6 1+ fe 3 )2-fe2 (b 2 +fc 3 )2-fe2 



i— — <5(w - b t - b 2 - k 3 



+ 



k ± k 2 k 3 (fc 1 +fc 3 ) 2 b\ {k 2 + k 3 ) 2 b 2 



5{uj — k 1 — k 2 — fc 3 ) 



(E19) 



The last term inside the large braces integrates to a result which is finite for arbitrarily small uj. The first six terms, 
which are multiplied by 5(lj — b 1 — fc-J, are finite for uj > m. The seventh term has a threshold uj = 2m. 



c. The O(g) contribution ac 
Proceeding in a similar manner, we obtain 



1 1 1 



+ 



b^b, (b 1+ b 2 )2-bl (b 1+ k 1+ b 2 Y-kl b^b, (fc 1+ b 3 )2-fc 2 2 (b 1 -b 3 )2-62 



<5(w — 6-l — k x ) 



k,b 2 b 3 (fc 1+ b 3 )2-fe2 (b2+b3 )2_ b 2 



5{uj — k 1 — b 2 — b 3 



(E20) 



The first Dirac delta has a threshold at uj = m, and the second Dirac delta has a threshold at u> = 2m. 

I 



where 



x(fc 1; k 2l ui) = 1 + 



w 2 — 2(fc 1 + fc 2 ) w 



2 /c-^ 



(E24) 



3. subgap conductivity to order g 

The complicated expressions derived above can, in 

principle, be evaluated numerically. Here we focus on the 

limit u < m, where we can obtain an analytic expression The constraints over the fc and k integrals are the fol- 

for the subgap conductivity as a power series in u/m. Of , , . 7,7,1 - 1 

lowing, .birst, since ui = k 1 + k 2 + k 3 , we must have 

k\ + k 2 < uj. Second, we must have x < 1, which gives 
k l -\-k 2 > hui. Finally, we must have x > — 1, which gives 
{2k l — ui)(2k 2 — uj) > 0. Putting this all together, we 
find that k 1 and k 2 are to be integrated over the shaded 
triangle in Fig. [9j 



all the terms computed above, only the first term in er A 
and the last term in a B are nonzero for ui < m. These 
add up to 



o- fl (w) 



Tig I d k l I d k 2 5((jj — k 1 — k 2 —k 3 ) 



4Nd J (2-n) d J (2Tf) d uik 1 k 2 k 3 m 2 +k 2 ~(k 2 +k 3 ) 2 



(N-l) kl k\+k\-k% 



2 +k 2 -(k 2 + k 3 ) 2 m 2 + k 2 -{k 1 +k 3 ) 2 



It is convenient to define k l = ^uiu and k 2 



^UJV. 



(E21) Then the vertices of the triangle in the (u, v) plane are 
(1, 0), (1, 1), and (0, 1). We then have 



We can write 



A/<^ — \l k-^ ~\~ k<2 ~\~ ^xk-^ k<2 •) 



(E22) 



-i r 

-(k.,+k,) 2 L 



kj(N-l) 



where x is the cosine of the angle between k 1 and k 2 . 
Then 



m 2 +k 2 -(k 2 +k 3 ) 2 l m 2 +k 2 -{k 2 +k 3 ) 2 m 2 + k 2 -(k 1 +k 3 ) 2 



uj 2 /4 r « 2 (jV-l) . 2(tiu-2ti-2t>+2) ] frpnr^ 
m 2 -w 2 (l-«) Lm 2 -w 2 (l-«) + m 2 -w 2 (l-t)) J l^^J 



<5(w - k 1 — k 2 — fc 3 ) 8{x — x(k 1 ,k 2l uj)) 



as well as 



U)k\k\ 



(E23) 



sc = H (1 — u — v) 



(E26) 
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FIG. 9. The range of (fci, fe) integration is the gray triangle. 

We obtain, in d = 3 space dimensions, the low- 
frequency behavior (ui <C m) 

u 2 (N-l) 



CT S M = 2^W^ du dv m 2 -J(l- M ) 



qe m 
3?r 2 2 9 



m 2 — cj 2 (1 — u) 

2(mj-2it-2iH-2) 
m 2 — cu 2 (1 — v) 

(E27) 

| u a *\ _,_ 9iV-16 ui 7 , 1 
N 180m 7 "T • • • J 



4m 3 10m- 



For d = 2 we find 



°"»( w ) " 2W Jdu Jdv y (u+ „_ 1) 1 (1 _ u)(1 _„ ) m »_J(i_„) : 



u J (JV-l) . 2(md-2m-2»;+2) 
i 2 — cj 2 (1 — u) m 2 — u 2 (l — i)) 



2 8 7T 



f AT-2 / 1 
I W I 1 



16w , 32w J \ , 3N-5 16. 



5m t 105m 3 j ~ N 315m 5 



Remarkably, in the case TV = 2, the conductivity van- 
ishes for small powers of w/m for both d = 2 and <i = 3, 
leading to a strong pseudogap behavior in the optical 
conductivity. 



Appendix F: Layered lattice bosons with Coulomb 
interactions 

The Lagrangian density for layered charged bosons is 
[see RefP, Eqs. (7-9)] 



Ant ' 



(Fl) 



where Q 2 = 16ire 2 /e is the effective electric charge of 
the bosons and C — diag(l, 1, a), with a 2 <C 1 the ra- 



tio of c-axis and in-plane superfluid stiffness. The bare 
propagators are now 

9 



G a a (fe) 



f P (k) 



< + Kb + « 2fc c + 1 



gQ 2 



fc2( w 2 +w 2 / 2 (A;)) 

l k 2 ab + a 2 p k 2 

V k 2 



(F2) 
(F3) 

(F4) 



Here oj p and a p u p are the in-plane and c-axis plasma 
frequencies given by 



Q 2 



<ib 



(F5) 



where are the lattice constants. 

We evaluate the O(g ) (one-loop) conductivity: 

x G^{iui m + iv n , kab, fc c ) (F6) 
We set [i = v = x and q = (iu) n , q = 0). We then have 



-(>■ 



Ef d 2 k ab f dkc fcq t 
J (2tt) 2 /c 2 J 2ir k l b +<xkl+l-(iv„ 



_Q 

2/3 



(F7) 



c fe2 {iv m +a k )(iv m -a k ) 

(F8) 



(ii> m -|-iu; n ) +iu n — b k ) 



where 



kl b + akl + l, b k =u p f p {k) (F9) 



(E28) This is of similar form to Eq. (E2), but with a different 



choice of a k and b k . We can then use Eqs. (E4) and (E5) 
to obtain the conductivity for ui > 0, 

CT M = & J dk *b J dk -^^ S - «fc - b k) (F10) 
Define k a b = kcosO, and anisotropy functions 
f(6) = Vcos 2 6 + a 2 sin 2 9 
fviO) = 

such that 



a 2 sin 2 i 



(FIX) 



a k = y/l + k*p(0), b k = co p f p (9) (F12) 

The conductivity integrals can be evaluated in cylin- 
drical coordinates, 



OO /' 7r /2 

dkk 2 / rf6> . cos36> . x 

x*(w- Vl + fc 2 / 2 (0)-w p / p (0)) (F13) 
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Changing the argument of the delta function from to to Thus, 
k 2 , provided (w — uj p f p {9)) > 1, yields 



(uj-u p f p (0)) 2 -1 
P(0) 



dk 2 

6(oj - c fc2 ) = 5{k 2 - fcS)"^f (F15) 



This yields a simplified expression, 



x e(w-w p /p(0)-l) (F17) 



This expression is plotted in Fig. [5] for the parameters 
(F14) o; p = 10, uj c = 0.1, and a = 10~ 2 . 
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